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LARGE-TIME BEHAVIOR OF PERIODIC ENTROPY 
SOLUTIONS TO ANISOTROPIC DEGENERATE 
PARABOLIC-HYPERBOLIC EQUATIONS 

GUI-QIANG CHEN AND BENOIT PERTHAME 

Abstract. We are interested in the large-time behavior of periodic en- 
tropy solutions in L°° to anisotropic degenerate parabolic-hyperbolic 
equations of second-order. Unlike the pure hyperbolic case, the nonlin- 
ear equation is no longer self-similar invariant and the diffusion term 
in the equation significantly affects the large-time behavior of solutions; 
thus the approach developed earlier based on the self-similar scaling 
does not directly apply. In this paper, we develop another approach 
for establishing the decay of periodic solutions for anisotropic degener- 
ate parabolic-hyperbolic equations. The proof is based on the kinetic 
formulation of entropy solutions. It involves time translations and a 
monotonicity-in-time property of entropy solutions, and employs the 
advantages of the precise kinetic equation for the solutions in order to 
recognize the role of nonlinearity-diffusivity of the equation. 



1. Introduction 

We study the large-time behavior of periodic solutions in to nonlin- 
ear anisotropic degenerate parabolic-hyperbolic equations of second-order. 
Consider the Cauchy problem for the second-order equations: 

(1.1) dtu + V^-f{u) = V^-{A{u)V^u), xGR^ t>0, 

(1.2) u\t=o = uo€ L°°{R'^), uo{x + P) = uo{x) a.e.xGR'^, 
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where P = {Pi, Pd) is the period with Pi > 0, Tp := Uf^^[0, Pi], / : M 
U.'^ satisfies 



and the dx d matrix A{u) = (ajj(u)) is symmetric, nonnegative, and locahy 
bounded so that it can be always written under the form 

d 



with {aik{u)) the square root matrix of A{u). 

Equation (II. Ih and its variants model degenerate diffusion-convection mo- 
tions of ideal fluids and arise in a wide variety of important applications (cf. 
[U [21 m [12 [20] and the references cited therein) , for which a deep under- 
standing of solutions to (jl.ip is in great demand, at both the theoretical and 
numerical level. 

In [7j, a well-posedness theory has been established for solutions of 
the Cauchy problem (jl.ip and (|1.2p of anisotropic degenerate parabolic- 
hyperbolic equations of second-order. It extends the isotropic theory for 
degenerate parabolic-hyperbolic equations and several latter studies, see for 
instance [3l [5l [TBI [T8] and the references therein. A notion of kinetic solu- 
tions, a new concept in this context, and a corresponding kinetic formula- 
tion have been extended. In particular, it has been also proved that, when 
uq G L°°, the kinetic solution is equivalent to the entropy solution, or to the 
dissipative solution [23], which is unique; this provides a new path to study 
the behavior of entropy solutions in through the corresponding kinetic 
equations for anisotropic degenerate parabolic-hyperbolic equations. In this 
paper, we employ the advantages of this path to develop a new approach 
for establishing the decay of periodic entropy solutions in to (jl.ip and 
()1.2p as t — > oo. The main theorem of this paper is the following. 

Theorem 1.1 (Main theorem). Assume that (jl.3p and (jl.4p hold. Let 
u G -L°°([0,oo) X W^) he the unique periodic entropy solution to (jl.ip and 
(jl.2p . Assume that the flux function and the diffusion matrix A{u) satisfy 
the nonlinear ity-diffusivity condition: For any 5 > {), 



(1.3) 



«(•) :=/'(•) Giioc(K;K') 



(1.4) 





k=l 



(1.5) 




\di 



UJS{\) 



0. 



sup / 

|r| + |K|>(5 J\£,\<\\u^ 



,0 



oo 



\ + \T + a{i)-n\^ + {K^A{i)KY ■ 



Then we have 



(1.6) 



t— >oo 







where 
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The nonlinearity-diffusivity condition (jl.Sp for equation (jl.ip is developed 
from [13^ [TT] and is reminiscent from the theory of velocity averaging lemmas 
[m EU [22l [25] . For smooth coefficients, condition (|1.5p is equivalent to the 
simpler and more standard setting: For any (r, k) G M*^"^^ with + |kP = 1, 
we have 

G M : lel < Ikolloo, T + a(e) ■K = 0, k^A{^)k = O} = 0. 

Several explicit examples are given in [17^ I25j . It implies that there is no 
interval of ^ on which both the flux function /(^) is affine and the diffusion 
matrix A{(^) is degenerate, and thus also makes the relation with the applica- 
tions of the theory of compensated compactness [191 [26] to one-dimensional 
hyperbolic conservation laws. Such a nonlinearity-diffusivity condition is 
necessary for the decay of periodic solutions and the compactness of solu- 
tion operators. See also the books [51 [2l]. 

Theorem 1.1 also extends naturally the non-degeneracy condition for the 
purely hyperbolic case in [11] where the first long-time convergence result 
of periodic solutions was obtained in one or two dimensions for BV initial 
data and higher order local non-degeneracy conditions that replace (jl.Sp . 

Unlike the pure hyperbolic case, equation (jl.ip is no longer self-similar 
invariant and the diffusion term in the equation significantly affects the large- 
time behavior of solutions; thus the approach in [6] based on the self-similar 
scaling for the pure hyperbolic case does not apply and we have to change 
the strategy of proof. The approach developed in this paper is based on the 
kinetic formulation of entropy solutions in [7] , involves time translations and 
a monotonicity-in-time of entropy solutions, and employs the advantages of 
the kinetic equations of the solutions, in order to recognize the role of the 
nonlinearity-diffusivity of the equation. 

The rest of this paper is organized as follows. We first recall the notion 
of entropy solutions and their precise kinetic formulation, and then ana- 
lyze some basic properties of entropy solutions. Finally, we develop a new 
approach to give a rigorous proof for the long-time asymptotic result. 



2. Entropy Solutions and Kinetic Formulation 

In this section, we first recall the notion of entropy solutions and their 
precise kinetic formulation, which requires some care to define appropriately 
the various terms of the equation. Then we analyze some basic properties 
of entropy solutions, which will be used in the proof of Theorem ll.il 

Definition 2.1. An entropy solution is a function u{t, x) G L°° ([0, oo) x R'^) 
such that 

(i) Edcc,Pikiu)£L^[0,oo)xR'^),k = !,■■■ , d, for Afc(^z) = paik{v)dv; 
1=1 
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(ii) For any function ^p G Co(M) with ip{u) > and any k = 1, - ■ ■ ,d, the 
chain rule holds: 

(2.1) 

d d 

Y.d^^l3l{u) = ^W)Y.^.Ak{u)^L\%^)ycM.''), for {(ifj = ^I3[^- 

i=l i=l 

(iii) For any smooth function S{u), there exists an entropy dissipation mea- 
sure m^"{t,x) satisfying that 

Jr 

such that 

(2.3) dtS{u) + V, • /(n) - V, • {A{C)V,S{u)) = -{m^" + n^") 

in P'(M+ X R'^) with initial data S{u)\t=o = S{uo), where g'^ : M ^ M*^ is 
the corresponding entropy flux, i.e. {q^{u)y = S'{u)f'{u), and 

n^"(t,x,e)= J S"{Onit,x,C)dC 

for n{t, X, ^) the parabolic defect measure of u{t, x) defined as: 

d d 

(2.4) n(t,x,e) :=<^(e-^^(i,x)) J^(^9,,Afc(n(t,x)))' a.e. 

k=l 1=1 

We point out that the chain rule in (ii) has to be assumed in the anisotropic 
case, and this makes the main difference with the isotropic case in |3] where 
this property follows from an argument reminiscent to the theory of Sobolev 
spaces. The requirement (iii) has been made with notations that are adapted 
to the kinetic formulation we introduce now. The L^-condition in (i) is re- 
quired to define the parabolic defect measure n(i,x,^) in (|2.4p (also see 
0). 

To do so, we may factor out an S'{u) in equation (j2.3p and obtain a more 
handful kinetic formulation of nonlinear degenerate parabolic-hyperbolic 
equations of second-order with form (jl.ip . The new ingredient of this for- 
mulation is the identification of the kinetic defect measure m{t, x, ^) and the 
degenerate parabolic defect measure n{t,x,^) in a precise manner, even in 
the region where u{t, x) is discontinuous. Compare with the classical kinetic 
formulation for scalar hyperbolic conservation laws in [17] (see also |21j). 

We introduce the kinetic function x on M?: 

( +1 for 0<^<u, 

(2.5) X{^;u) = < -1 for n<^<0, 

I otherwise. 

We notice that, if u is in L°°([0, oo) x M'^) and periodic in x with period P, 
then x{S,', u) is in L°°([0, oo) x M'^+^) and periodic in x with period P. 
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The simple representation formula S{u) = S' {^) xiC', u) leads to the 
following kinetic equation, which is equivalent to the entropy identity (j2.3p : 

(2.6) dtx{C;n) + a{0-^MC;u)-V^-{A{C)V^MC;u)) = d^{m + n){t,x,C) 
in V'{R+ X M'^+i) with initial data 

(2.7) xiC;u)\t=o = xiC;m), 

where n(t, x,^) is still defined through (12. 4p . 

In [7], it has been proved that the entropy solutions in L°° are equivalent 
to the kinetic solutions determined by the kinetic formulation (|2.5p - (j2.7p . 
Furthermore, we have 

Theorem 2.1. Assume that (fOP and ([Ll]) hold. Then 

(i) there exists a unique entropy solution u E -L°°([0, oo) x M*^) to (jl.ip ~ (|1.2p 
such that u{t, x + P) = u{t, x) a.e. 

(ii) the entropy solution u{t, x) satisfies 

(a) \\u{t, < ||no||L°o; 

(b) for any t2 > ti > 0, 

(2.8) / \u{t2,x)\^ dx < / \u{ti,x)\^ dx, 

JTp JTp 

(2.9) / \u{t2,x) — v\dx < / \u{ti,x) — v\dx for any constant v . 

JTp JTp 

(iii) Moreover, if the flux function f{u) and the diffusion matrix A[u) further 
satisfy the nonlinearity-diffusivity condition (II. 5|) . then the solution operator 
u{t, •) = S'jno(-) : L'^ — > is locally compact in {t,x) for t > 0. 

The results (i)-(ii) in Theorem 2.1 are direct corollaries of the well- 
posedness results and the arguments in [7] (i.e. standard entropy inequalities 
for |u — t;| and u'^) as in the hyperbolic case [8| [24]; while (iii) is a direct 
corollary of the kinetic averaging compactness result of [I7j; also see more 
recent results in [T5l[22l[25] . 

3. Decay of Periodic Entropy Solutions: Proof of Theorem 11.11 

In this section, we develop a new approach to give a rigorous proof for 
the decay property of periodic solutions, which takes the advantage of the 
precise kinetic formulation (12.5p -( f27fll . Without loss of generality, we first 
set Jj.^ u{t,x)dx = 0; otherwise, we may replace u{t,x) by u(t,x) — u, f{u) 
by f{u + u), and A{u) by A{u + u) in (jl.ip . so that all the arguments below 
remain unchanged. Then, we divide the proof into four steps. 

1. Limit. Theorem 2.1 indicates that the periodic solution u{t,x) belongs 
to L°°, bounded by ||no||Loo, and is compact as the solution operator. Also 
the function 

I{t) := / \u{t,x)\'^dx 

JTp 
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is a non-increasing, bounded function, which imphes that the following limit 
exists: 

(3.1) lim I{t) = /(oo) =: £ [0, oo). 

t—KX) 

2. Translations. Set 

Vk{t, x) := u{t + k,x). 

Then we find that, for t > —k, 

(i) \\vk{t,-)\\L°° = \\u(t + k,-)\\L°° < \\uo\\l°^; 

(ii) Vk{t,x) is also a periodic entropy solution to (|l.ip - (|1.2p : 

(iii) for each A; > 0, xi^jV^i^jx)) satisfies 
(3.2) 

dtxiC; v") + aiO • VMC; v") - V. • (^(OV.x(^; v')) = di^{m^ + n^){t, x, 
in V'{{-k,oo) X M'^+i). 

By Theorem 2.1 (iii) applied to v'^, there exists a subsequence {f'^^ j^i C 
and v{t,x) £ L°°(M'^+^), with J^^ v{t,x)dx = 0, such that 

v'^^ (t, x) v{t, x) a.e. (t, x) G as j oo. 

Correspondingly, we have 

x{C;v''{t,x)) ^ X{tv{t,x)) a.e. {t,x,^) as j ^ oo. 

Furthermore, multiplying both sides of (|3.2p by ^ and then integrating 
over (t,x,0 e {-T,T) x Tp x R for any T > 0, we obtain 
(3.3) 

/ / {m^ + n^){t,x,i)dtdxdi<]-{l{k-T)-I{k + T)) <\\Tp\\\u4l^. 
J-T JTp ^ ^ 

This implies that the nonnegative measure sequence (m^ +n^)(t, x,^) is 
uniformly bounded in k over {—T,T) x Tp x M, and hence there exists a 
subsequence kj and a measure M(t, x, ^) such that 

(m^^ + n^^)(t,x,^) ^ M(t,x,^) > weakly in as j ^ oo. 

On the other hand, since I{t) converges, we also have 

(3.4) I{k-T)-I{T + k)^0 as k ^ oo. 

Then we conclude from (j3.3p that M((— oo,oo)xTpxM) = 0, which implies 

(3.5) M(M'^+2) = 0. 

Furthermore, letting A; ^ oo in (|3.2p . we conclude that x(Cj ^) is a Tp- 
periodic solution in D'(M x 1^"^+^) of 

(3.6) dtx + aiO ■ V,x - • {A{OV,x) = 0. 

In particular, multiplying (j3.6p by and then integrating dxd^, we have 

(3.7) / |f(t,x)|^dx = /oo e [0,oo), VtGM, 
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where I^o = is a constant, independent of t, determined in (13.11) . 

The rest of the proof consists in showing that such a function x is very 
particular and is in fact constant (see also this type of "rigidity" results in 

3. Claim: v{t,x) = a.e. for {t,x) eRxR'^. 

We introduce a time truncation function (j){t),0 < (p{t) < 1, so that (px 
belongs to ^^(k x Tp x M). Then we have 

(3.8) dt{M+a{0-^A<Px)-^.-iA{0^,{cl)x)) =xdt^ mV'{R''+^). 

Next, since (px and x4't are periodic in x, we take the global Fourier 
transform in t G M and the local Fourier transform in x G Tp to obtain 
g{T,K;^) for ((/>x)(*; 2;, ^) and /i(r, k;^) for {xdt(j)){t,x,S,) in L^, where the 
frequencies k = {ki, • • • , k^) are discrete: 

27r 

i^i-B-, n = 0,±1,±2, • • • . 

-< i 

That is, for example, 

5(T,K;e) = ^/ / {M{t,x,Oe-'^-'+^-^Utdx 

I J^pI Jr JTp 

so that 

(cPx)it,x,0 = ^Y. I 5(r,^;e)e*(^*+'^-^)dr. 

Taking the global Fourier transform in t E M and the local Fourier trans- 
form in X G Tp in both sides of (|3.8|) . we obtain 



' i{T + a(0 • n) + K^A{S,)Kjg = h. 

Following usual ideas from the kinetic averaging lemmas, we may introduce 
a free parameter A > (to be chosen later on) and write 

(yX + i{T + a{(,) ■ k) + tJA{^)K^g = h + VXg. 

This leads to 

Q = (h + V^g) —F= : ^ . 

^ VX + i{T + a{0-K)+K^AiC)K 

Integrating in ^ and using the Cauchy-Schwarz inequality, we find 
<2( [ h^di + X I g^d^ 



2 

di. 



\/A + i(r + a{£) ■ k) + A{£,)k 

Notice that the frequencies k are discrete and may include k = 0. In 
particular, when k 7^ 0, then there exists (5o > such that \k\ > 5q. On the 
other hand, since v{t,x) has mean zero in x over Tp, we have 

^(r,0) = 0. 
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The nonlinear ity-diffusivity condition (jl.Sp gives that, when «; / 0, for 
any 6 G (0,(5o), 

\^\^<C'^ [ h^d^ + Cus{\) [ \g\^d^, 



and thus 





A 



<C^^4^ I ix^tfdtdxdC + CiOsiX) I \(l)x\'dtdxdi. 



This implies 

f \(l)v\'^dtdx<c'^^^ [ \(l)t\'^\v\dtdx + CuJsiX) [ \(l)\'^\v\dtdx. 

Using (j3.7p and the Cauchy- Schwartz inequality, we arrive at 

(3.9) loo fu\(l>\^dt <Cus{X)( [ Ivl^dxV^W [ \^t\^dt+ I 

^ Jtp ^■^ Jr Jr 



k\'^dt+ I i^i'^dt]. 



Choosing first A small and then small, we conclude from (|3.9p that 

/oo = 0, 

which implies from (j3.7p that 

(3.10) v{t, x) = a.e. {t, x) G M x R'^. 

On the contrary, if I^o > 0, then we can choose A small enough so that 
CiOsiX)/lU'^ < i and find from ^M) that 



[ \^\^dt < 2ClU^^^ [ l^tl^dtdxdC. 
Jr X 



It remains to choose a sequence of functions (j)Bit) = 1 for \t\ < B, with B 
a given large number and (p'si^) ~ ioi B < \t\ < 2B, and (psit) = 

for |t| > 2B. In the above inequality, we find 

t1/2 <■ ^ ^sjX) 

°° - 52A ' 

where C > is a constant independent of B and A. When B tends to oo, 
this implies that /qo must vanish, which is a contradiction. 
Therefore, (I3.10p holds. 
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4- Conclusion. Then 

\v''^{s,x)\ dxds — > as J — > oo. 

/o JTp 

Therefore, for any T > kj + 1 > kj + s iov s ^ (0, 1), we employ ()2.8p for the 
monotonicity-in-time of solution to obtain 

'1 . ^ 



'0 JTp 

We conclude that 



v'''{t,x)\dxdt> I \u{l + kj,x)\dx> I \u{T,x)\dx. 



/ \u{T,x)\dx^^ as T ^ oo. 
This completes the proof. 
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